LIBS-Acoustic Mid-Level Fusion Scheme for Mineral Differentiation under Terrestrial and Martian Atmospheric Conditions

The shockwave produced alongside the plasma during a laser-induced breakdown spectroscopy event can be recorded as an acoustic pressure wave to obtain information related to the physical traits of the inspected sample. In the present work, a mid-level fusion approach is developed using simultaneously recorded laser-induced breakdown spectroscopy (LIBS) and acoustic data to enhance the discrimination capabilities of different iron-based and calcium-based mineral phases, which exhibit nearly identical spectral features. To do so, the mid-level data fusion approach is applied concatenating the principal components analysis (PCA)-LIBS score values with the acoustic wave peak-to-peak amplitude and with the intraposition signal change, represented as the slope of the acoustic signal amplitude with respect to the laser shot. The discrimination hit rate of the mineral phases is obtained using linear discriminant analysis. Owing to the increasing interest for in situ applications of LIBS + acoustics information, samples are inspected in a remote experimental configuration and under two different atmospheric traits, Earth and Mars-like conditions, to validate the approach. Particularities conditioning the response of both strategies under each atmosphere are discussed to provide insight to better exploit the complex phenomena resulting in the collected signals. Results reported herein demonstrate for the first time that the characteristic sample input in the laser-produced acoustic wave can be used for the creation of a statistical descriptor to synergistically improve the capabilities of LIBS of differentiation of rocks.

L aser−matter interaction encompasses a large number of phenomena taking place over a time lapse from femtoseconds to milliseconds. 1,2 The rapid succession in which events unfold after the laser light reaches the sample turns this process into a complex scientific problem. Still, as demonstrated by decades of research, thoroughly comprehending the different stages from the absorption of laser photons and changes in the optical traits of surface to thermal relaxation, material ablation, and plasma formation is of great interest for analytical applications. Laser-induced breakdown spectroscopy (LIBS) takes advantage of the elemental and molecular emission to reveal physicochemical information pertaining to the inspected sample. 3,4 Currently, LIBS is a mature technique that can face the chemical characterization of a variety of samples inaccessible to a majority of other analytical tools. 5 To name a few figures of merits, the scarce amount of sample required allowed the analysis of single nanoparticles with masses down to the attogram scale. 6 The capability to carry out simultaneous multielemental characterization with a high spatial resolution has been exploited for a wide range of targets including biological, industrial, and geological samples. 7 Additionally, fieldable and remote sensors are a key application niche for LIBS where it truly exhibits its adaptability, 8 from submerged shipwrecks, 9 through geochemical in-field analysis, 10 to space exploration. 11 Shockwave generation is one of the phenomena occurring alongside the laser-induced plasma owing to the explosive ejection of ablated material into the surroundings. The plasma shockwave is expected to encompass two different, entwined contributions: one related to the expansion of the plasma, and another linked to the mechanical relaxation of the probed material. Despite the interest shown previously in exploring laser-induced acoustics as a means to expand the information delivered by LIBS inspection of solids, the multiple sources of uncertainty linked to sound waves have hindered their potential application in sensors for both in-lab and off-lab studies. Nonetheless, the recent integration of a microphone synchronized to the LIBS laser in the SuperCam instrument in the Perseverance rover deployed as part of the NASA Mars 2020 mission has re-kindled the motivation toward unraveling the shockwave as it evolves in time to extract the samples' contribution from it seeking to ascertain parameters such as the crystalline structure and the presence of layers or to enhance the overall SuperCam performance. 12−14 Early reports on laser-induced acoustics outlined the relation existing between acoustic energy and the energy density at the sample 15 and the possibility to discriminate the mechanism leading to plasma formation on the basis of the recordings with further research being focused on using acoustics to normalize LIBS data. 16 More recently, a series of studies by Murdoch et al. and Chide et al. set the basis for modern-day laser-induced acoustics while performing tests for the SuperCam instrument. 14,17 These authors carried out the recording of plasmas formed under Earth and Martian atmospheric conditions (i.e., about 7 mbar CO 2 , T down to −80°C at the microphone and different intensity wind currents) on Martian simulant targets as well as in minerals of interest for the mission. Results indicated a strong dependence of the generated acoustic wave with physical properties such as the hardness of the material, and therefore the ablation rate, the density, and the thermal conductivity. Sound propagation as a function of the medium and the length of the acoustic path was also evaluated in refs 17, 18 since the sound wave is known to strongly attenuate with distance and the LIBS + acoustic tandem is intended to be used to measure samples located 2−7 m away from the rover. In spite of positive results indicating the complementarity of acoustics and additional information it can provide, many questions still remain as to how LIBS data, specifically for the identification and classification of mineral phases, can be improved by merging of simultaneously acquired sets.
In the present work, Fe-based and Ca-based minerals featuring intraseries similar LIBS spectra are probed at a sampling distance of 2 m under Earth and Mars-like atmospheres. LIBS data and simultaneously acquired acoustic recordings were pre-processed to develop a mid-level data fusion strategy, based on a combination of LIBS principal components analysis (PCA) and acoustic features to create a new sample descriptor allowing better differentiation of samples with extremely similar LIBS spectra. After the model was verified for mineral data acquired in Earth atmosphere, the training set was sampled under Mars-like conditions. Results indicate that, under a tightly controlled experimental scenario, acoustic data in the time domain can be merged with LIBS spectra to enhance the discrimination capabilities of the technique owing to differences found in the laser-produced sound wave, thus paving the way toward extracting relevant chemical data from this newly exploited source of information.
■ EXPERIMENTAL SECTION Experimental Setup. Experiments reported herein were performed in both terrestrial and Martian atmospheres. Measurements were conducted inside a thermal vacuum chamber (TVC) installed at the UMALASERLAB ( Figure  1a). This facility was built by AVSAdded Value Solutions (Gipuzkoa, Spain) for remote spectroscopic studies in Earth and extraterrestrial environments. The TVC is a 12 m long and 2 m diameter cylinder with a total useful volume of ca. 24 m 3 capable of recreating the atmospheric conditions of different planetary bodies within a range of 200−400 K and 10 −3 −10 3 mbar. 19 The chemical composition of the working atmosphere could also be regulated by a gas inlet using a mass flow controller.
A schematic diagram of the experimental setup is presented in Figure 1b. A pulsed Nd:YAG laser (Ekspla, model NL 303D/SH; @1064 nm, 45 mJ, pulse width 4 ns) located outside the TVC was used as the excitation source. Collimated laser pulses were expanded and then focused using a beam expander prior to entering the TVC through a borosilicate glass vacuum viewport. Once inside, laser pulses were guided toward the target (Figure 1c-1) for generating the recorded plasmas. Samples were placed upon a platform coupled to a motorized linear stage allowing displacement in the axis perpendicular to the laser path ( Figure 1d). Plasma emission was remotely collected using a 6′ Ritchey−Chretien telescope (f/9, focal distance 152 mm) modified to integrate an optical fiber (Figure 1c-1−3). The plasma image was coupled to a 1000 μm single optical fiber, which was bifurcated after the fiber optic vacuum feedthrough to SMA terminated 600 μm fibers, with a total length of 3 m. A pair of miniature Avantes Czerny-Turner spectrographs with a 75 mm focal length placed outside the TVC was also used with a fixed integration time of 1.28 μs. Under this configuration, LIBS spectra ranging from 240 to 590 nm were recorded.
Acoustic waves coming from laser ablation and plasma expansion were recorded using a 6 mm pre-polarized condenser microphone (20 Hz to 19 kHz frequency response, omnidirectional polar pattern, 14 mV·Pa −1 sensitivity, TR-40 model from Audix). The microphone (Figure 1c-4) was housed inside a custom-built hemi-anechoic chamber (50 cm length, 25 cm internal diameter) with inner and external walls insulated using acoustic foam. This acoustic suite was positioned at a fixed sample surface-to-microphone distance of 2 m. The source−receiver path featured an angle of about 7°t o avoid obstructions in the optical path. A 24-bit/192 kHz audio interface (UA-55 Quad-capture, Roland) was used at a sampling rate of 96 kHz for the digitalization of acoustic waves. Audacity software was employed as an audio recording application. The sample holder was placed inside another custom-built hemi-anechoic chamber (70 × 40 × 40 cm 3 , L × W × HFigure 1d) covered with HiLo-N40 acoustic foam made of polyurethane with high rigidity and low density (70 mm total thickness, 40 mm knob height, 16.5 kg·m −3 bulk density). This optical configuration led to a spot diameter of ca. 300 μm at the surface of samples located at 2 m distance away from the last mirror and the collection optics. To avoid acoustic reflections, the ground between the samples and the microphone was also covered with anechoic foam.
Samples. To evaluate the input of the acoustic response to the information yielded by laser-induced plasmas and better identify geological specimens, 12 minerals, 6 rich in Fe and 6 rich in Ca, were selected. These minerals were chosen as their expected elemental composition should translate into very similar LIBS spectra, either due to the presence of some components in the atmosphere, thus making its source ambiguous, or due to their difficult detectability by LIBS (e.g., S).
Pyrite (FeS 2 ), siderite (FeCO 3 ), magnetite (Fe 3 O 4 ), hematite (Fe 2 O 3 ), limonite (FeO(OH) n H 2 O), and goethite (FeO(OH)) were selected as Fe-bearing minerals. Iron is among the 18 most abundant elements in our solar system, and, particularly considering the interest in Mars exploration, these specimens have been both detected in Martian meteorites and identified in situ in several regions on the Red Planet. 20−22 Moreover, a group of six calcium-rich minerals was evaluated. Calcium is a key constituent of rock-forming minerals and can yield clues about the material that formed the Solar System's rocky planets (Mercury, Venus, Earth, and Mars). In this context, limestone, a sedimentary rock consisting primordially of calcium carbonate, alongside aragonite, calcite spar, and calcite, all with the same expected composition (CaCO 3 ), were probed owing to their close LIBS response. Finally, gypsum, common sulfate mineral composed of hydrated calcium sulfate (CaSO 4 ·2H 2 O) and its fine-grained massive variety, called alabaster (CaSO 4 ·2H 2 O), were also investigated. 22−26 Prior to their analysis, the minerals, originally in their natural form, were cut and polished. A total of 75 laser events were analyzed, i.e., 25 plasma events per position at three different positions along the surface of the mineral.
Data Processing. A new global descriptor was created to perform data fusion via a combination of the spectral responses from both the laser-induced plasma emission and the acoustic wave, to describe the interrogated targets. This combination can be done at different levels: low-levelworking with the spectra by simply concatenating or applying the outer vector over the responses; mid-level or intermediate fusionfor which some relevant features from each source are treated independently and then concatenated into a single array used Analytical Chemistry pubs.acs.org/ac Article for multivariate classification; and high-level fusionthe classification or regression models are run from each data array, and the results are combined to obtain the final result. 27 The second approach was applied in this work. It should be stressed that, although LIBS and acoustic responses are originated from the same laser event, the mechanisms governing each wave differ. In LIBS, the signal mainly derives from atoms (neutrals or ions) that have sequentially undergone a process of ablation, fragmentation, atomization, ionization, and excitation, in a proper temporal range, and at a wavelength range detectable by the spectrometer. In the case of acoustics, the main wave is attributed to the compression/rarefaction generated due to plasma expansion. For this reason, the combination of both data can provide complementary information to obtain this new attribute aimed to identify the interrogated specimen more clearly. Figure 2 shows the full data fusion scheme applied. In the intermediate data fusion scheme followed herein, different features were extracted independently from each emission and sound wave spectra data set. In the case of LIBS spectra, data were averaged to convert a set of 12 × 3 × 25 (samples × positions × spectra number) into a 12 × 3 one,  Analytical Chemistry pubs.acs.org/ac Article generating an average spectrum for each position. Therefore, an n × j matrix was obtained, where n corresponds to the number of median spectra and j corresponds to the number of spectral variables. These median spectra were normalized using unit vector normalization (eq 1) and mean centering (eq 2).
( ) Principal component analyses were performed using the resulting matrices. PCA provided the weights needed to get the new variable (principal component) 28,29 that best explained certain underlying facts from a multivariate data set without diluting essential information. Therefore, the PCA served not only as a variable reduction step but also as exploratory data analysis (EDA), being a zero step prior to the supervised classification model. 30 In the acoustic dataset, two different features were extracted from the time domain responses: the wave peak-to-peak amplitude and the variation of this value for a specific position. Both values are linked to the physical traits of the material, 14,17 supplying complementary information to LIBS spectra. Then, the PCA score matrix and the acoustic features were concatenated, obtaining the new descriptor in the mid-level fusion process. This descriptor was used as an input matrix for studying the classification's throughput for the minerals using linear discriminant analysis (LDA). LDA creates a linear boundary between classes based on the distance to each class' centroid. This is calculated using the Mahalanobis distance. 30 All of the data analysis described in this subsection was performed using Matlab 2020b, Mathworks, Inc. software.

■ RESULTS AND DISCUSSION
Exploratory Data Analysis on LIBS Spectra. First, exploratory data analysis on the emission spectra obtained by LIBS was performed. Figure 3 shows the average spectra for all of the samples (Figure 3a,d). As expected, LIBS information of minerals containing line-rich elements (e.g., Fe or Ca) exhibit almost identical spectral features. The minimal differences detected between collected LIBS data may be attributed to fluctuations in the morphological and physical characteristics of the material guiding the ablation rate and altering the lasermatter coupling rather than to changes in the concentration of the constituent elements. Moreover, as observed in Figure 3, the emission of some minor or unexpected elements was present in the spectra (Figure 3b,c,e,f). For Fe-rich minerals, these lines were mainly attributed to Mg, Ca, and Si, whereas for Ca-rich minerals, Mg and Sr were the most relevant impurities. As discussed in the Experimental Section, PCA was applied on the median LIBS spectrum for each position after continuum background subtraction, normalization, and mean centering. Figure 4 presents the score values for the first three principal components, with a total explained variance of 97%. From Figure 4, it is readily apparent how the two subgroups studied were clearly separated through the PCA score threedimensional (3D) graph. This was expected due to the notorious compositional differences found between the two subgroups. A partial, yet incomplete, separation of each mineral phase within each of the subgroups can be observed. This behavior is justified by the presence of impurities in each mineral, as shown in Figure 3, that create the cluster pattern observed in Figure 4.
Exploratory Data Analysis on Acoustic Data. The acoustic waves produced by the laser-induced plasma detonation were recorded simultaneously with the LIBS spectra. Figure 5a represents average spectra in the time Since an open-air path connected the two hemi-anechoic chambers, some of these reflections could be recorded by the microphone. Still, the different pathlengths traversed by reflected sound waves, echoes, and the original plasma blast were long enough to separate these contributions in time, thus allowing us to clearly temporally differentiate them and to choose only the portion directly related to the sample. Therefore, the information provided by the first maximum (t = 0 in Figure 5) of the acoustic wave should be completely free from interferences. As mentioned, two different values of interest were extracted from the peak-to-peak amplitude, its total value, and the coefficient of intensity variation as the material was drilled by the laser in one sampling position. This coefficient, given by the slope of the amplitude as a function of the number of shots, provides the change of the acoustic signal when drilling down into the sample due to successive laser shots. Figure 5b,c shows the variation of the average amplitude in the first maximum of the acoustic wave for each mineral studied ( Figure 5b) and of the slope (Figure 5c), showing a clear difference between the two subgroups. The values obtained in Figure 5 point toward a clear difference between the Ca and Fe groups, with similar values inside each group. Furthermore, the negative slope present for the Ca group indicated lower hardness and higher ablation rate per sample. 14,17 Data Fusion Treatment. The mid-level fusion step was carried out taking into consideration the descriptors obtained during the exploratory analysis described in the previous sections. From the LIBS point of view, four PCs were selected (>97% of cumulative percentage of the total explained variance). Therefore, the new fused descriptor was constructed by concatenating the score arrays from the LIBS-PCA with the values coming from the amplitude of the acoustic wave and the slope of its decay at each sampling position. This new descriptor is graphically represented in Figure 6, allowing us to clearly differentiate between the two main subgroups. Each group presented a characteristic graphical trait, with the Fe-rich minerals being the more intense in the acoustic amplitude and in the PC1 score. This PC showed more relevant presence of Fe lines in the PCA loadings, while the Ca-related scores had negative scores. PC2 and PC4 seemed to have similar values for all of the samples, exhibiting less subgroup-dependent behavior and, probably, reporting information for the intragroup sample clustering. The acoustic input here is equivalent to that seen in Figure 5a,b.
Once the new descriptor was built, it was evaluated to assess its feasibility to discriminate different minerals according to the fused optical and acoustic features. The classification performance of this new descriptor was, therefore, analyzed using an LDA model on the samples descriptors. The hit rate was individually calculated as follows for each data set: The data matrix, composed of the new descriptors built for the different sampling positions in each sample, was randomly divided into two matrices, i.e., the training matrix (75% of the samples) and the validation matrix (25%). Since there were only three sampling positions per sample, we imposed that there were no more than one replicate per sample in the validation matrix. By doing so, samples were not underrepresented in the training matrix. Subsequently, the LDA model was obtained using only the training matrix, and this model was then applied to the validation matrix to obtain the predicted group. These predictions were compared to the actual class values for the training matrix and the validation matrix. This process was iteratively repeated 1000 times with the resulting prediction values being the average of the percentage of correct assignations for all iterations. Due to the randomness of the matrix selection process, a limitation to prevent two different iterations from having the same components was also imposed. Table 1 shows the success rate for each data source, for the fused data and the original ones. LIBS success rate was obtained using four PCs (LIBS-PCA column) and the original LIBS spectra (normalized and mean-centered). Training and validation data are also represented, providing the training matrix better predictions, as expected. It is clear that fused data improved the classification capabilities upon comparison with the original data. Acoustic features on their own provided the worst hit rate between the studied cases. Still, the hit rate was  Analytical Chemistry pubs.acs.org/ac Article above 74% in all of the studied cases. Also, LIBS spectralrelated hit rate slightly overcame the LIBS-PCA, due to the presence of more qualitative information in the original spectra than in the PCA's score values. Furthermore, the intragroup success rate was also calculated, that is, using only samples belonging to one of the two inspected families being considered. Figure 7 shows the confusion matrix for the obtained results from the fused data. Figure S1 in the Supporting Information presents the independent confusion matrix for LIBS-PCA and acoustics data. For the Ca-rich samples, the more common false prediction was related to extremely similar mineral phases (i.e., alabaster-gypsum or aragonite-spar-calcite), while the main source of error in Ferich samples corresponded to the magnetite-limonite pair. Both minerals had similar average contributions from Mg, Si, and Ca elemental emission lines, so failures in the assignation may be attributable to these spectral presences. No false identification between the different subgroups was observed due to the dramatic differences in both emission spectral and acoustical responses reported. Application of the Data Fusion Strategy to Mineral Information Collected under a Martian Atmosphere. The preceding sections have shown that the mid-level fusion of optical emission and acoustic data from the laser-induced plasmas improves the classification capabilities for mineral phases in comparison to LIBS data alone. As a further step in this work, the developed strategy was applied to the same sample set under Mars-like atmosphere. Mars atmospheric pressure varies between 6 and 10 mbar depending on the geographical area and the season. 31 The composition of the Mars atmosphere was determined by the Viking landers in the 1970s. Carbon dioxide is the main constituent (95%), followed by nitrogen (2.7%), argon (1.6%), and others (<0.5%). 32 The parameters in the current experimental atmosphere were defined as 8 mbar of a CO 2 -rich atmosphere (>95%) with a sample temperature of 273 K.
The changes in the behavior of the laser-induced plasma caused by variable atmospheric pressure and composition have been widely reported; 33−35 however, the propagation of plasma-induced pressure waves is less known. Changes in the transmission of the acoustic wave (which will depend on P, T, and composition of the wave's path) were expected both in sound intensity and in the specific absorption of a range of frequencies.
Apart from the specified changes, other experimental parameters remained as described in Section 3. In this scenario, lower signals were detected for both LIBS spectra and the acoustic wave. Figure S2 in Supporting Information compares the LIBS emission spectra on Mars-like and Earth conditions for calcite and pyrite. Overall, less intense emission intensity was observed, as well as the signal-to-noise ratio (SNR) strongly decreased in the Mars atmosphere due to changes in the plasma dynamics, decrease in the electronic temperature, and thinner optical plasma density. Table 2 presents the SNR values for two selected emission lines, Ca II at 393 nm and Fe I at 388 nm, in calcite and pyrite, respectively. Moreover, it should be noted that the molecular emission bands were absent in calcite. The most prominent molecular systems featured in Earth LIBS spectra for calcite were CaO/CaOH at 550 nm. Apart from the mentioned physical changes, the pressure and the composition changed; thus, the drop may indicate recombination with dissociated atmospheric O 2 as the main source of the CaO band. This statement is further supported by the fact that the CO 2 found  Analytical Chemistry pubs.acs.org/ac Article in the Martian atmosphere has a higher dissociation energy than O 2 (532 vs 498 kJ mol −1 at 298 K), relinquishing a lower number of O atoms into the plasma. Following the procedure described above, a PCA was carried out on the LIBS data. The spectral matrix was also mean-centered, and the explained variance was obtained, yielding 95% for the first three PCs. Similarly, the Fe-rich and Ca-rich samples were grouped in the 3D scores plot (Figure 8).
Acoustic waves were also measured and processed as previously, exhibiting a considerably reduced signal intensity. Due to the pronounced decrease of the acoustic signal (almost 2 orders of magnitude) and the unavoidable decrease of the SNR (Table 2) in the simulated Martian atmosphere, most of the Ca samples, which had the lowest emission intensity in the previous study, did not provide a statistically detectable signal. Although this reduced the sample set, the mid-level data fusion  Analytical Chemistry pubs.acs.org/ac Article strategy allowed no cross-predictions in the class assignation between groups (Figure 9), so the study of the Fe-rich subgroup was expected to provide an idea of the behavior of the data fusion under different conditions. Limestone and calcite data were preserved, as their signal was larger than the microphone threshold. To improve the identification of the peaks by the algorithm, a bandpass filter was applied between 500 and 20 000 Hz, which eliminated part of the detected noise and allowed the identification of the signal coming from the plasma. Following the workflow presented in Figure 2, LIBS-PCA scores values and the acoustic features were concatenated and an LDA was applied in parallel to the procedure followed in the previous section. The success rate for the data fusion treatment was lower in Mars-like atmosphere than under Earth conditions (Table 3); however, the results obtained still showed a promising perspective for the fusion treatment. Lower SNR values were identified as one of the main causes of the lower hit rate given by the model besides the fading of some spectral features (such as the molecular emission) explained above. The confusion matrix for Martian data is presented in Figure 9. Following results in the terrestrial atmosphere, the cross-predictions between limestone and calcite were expected due to the composition of the samples. This same trend can also be found for the limonite− magnetite pair, as previously seen. Cross-prediction was also found for hematite-magnetite. In the work by Chide et al., 36 laser-induced phase transition from hematite to magnetite was observed; these thin layers along SNR collapse could lead to the generation of false predictions.

■ CONCLUSIONS
Herein, a data fusion strategy merging simultaneously acquired LIBS and laser-induced acoustic data was developed to improve the discrimination of spectrally similar minerals in a remote LIBS configuration. A new descriptor used in a midlevel fusion approach was able to improve the discrimination success rate, via LDA, compared to LIBS alone on the two mineral families studied, i.e., Fe-based and Ca-based rocks. The combination of LIBS-PCA and acoustics provided an easily explainable descriptor with a few variables, which avoided the handling of complex data arrays. This combination was based on the fusion of scores values obtained from LIBS data, and the identified acoustic features in time domain, viz., peak-topeak amplitude, and intraposition amplitude slope. Once validated under Earth conditions, the approach was tested for robustness under Martian conditions. Upon comparison to only LIBS or only acoustic datasets, the new descriptor also improved the discrimination hit rate in the Mars-like atmosphere. The success classification is increased from 90% (LIBS-PCA) and from 77% (acoustic) to 92% using the fused descriptor under Earth atmospheric conditions, while in Marslike atmosphere, the increment goes from 85/81% (LIBS-PCA/acoustics) to 89% (mid-level fusion). While this paper constitutes a preliminary attempt at fusing LIBS and acoustics to maximize the information yielded by analysis events, its possible implementation in the analysis of samples in an open environment must be carried out carefully, especially when acoustic data are involved. The presence of echoes or interferences can modify the time domain signal if the path difference is small enough to interfere (either constructively or destructively) in the prompt acoustic wave, thus changing the peak-to-peak values. One remaining issue is the slight decrease of the hit rate in Martian atmosphere data. Instruments with a better signal-to-noise ratio (more sensitive microphones or the use of an intensified CCD, for instance) could improve these results and bring them on a par with those obtained in the terrestrial atmosphere.
Confusion matrices for LIBS-PCA and acoustic inputs ( Figure S1) and comparative mineral LIBS spectra in Earth and Mars-like atmosphere ( Figure S2)    Analytical Chemistry pubs.acs.org/ac Article